Mie scattering by a charged dielectric particle 



R. L. Heinisch, F. X. Bronold, and H. Fehske 
Institut fur Physik, Ernst-Moritz-Arndt-Universitat Greifswald, 17489 Greifswald, Germany 

(Dated: November 26, 2012) 

We study for a dielectric particle the effect of surplus electrons on the anomalous scattering of light 
arising from the transverse optical phonon resonance in the particle's dielectric function. Excess 
electrons affect the polarizability of the particle by their phonon-limited conductivity, either in a 
surface layer (negative electron affinity) or the conduction band (positive electron affinity) . We show 
that surplus electrons shift an extinction resonance in the infrared. This offers an optical way to 
measure the charge of the particle and to use it in a plasma as a minimally invasive electric probe. 

PACS numbers: 42.25.Bs, 42.25.Fx, 73.20.-r, 73.25.+i 



The scattering of light by a spherical particle is a fun- 
damental problem of electromagnetic theory. Solved by 
Mie in 1908 [l| , it encompasses a wealth of scattering phe- 
nomena owing to the complicated mathematical form of 
the scattering coefficients and the variety of the under- 
lying material-specific dielectric constants 0, Q • While 
Mie scattering is routinely used as a particle size diagnos- 
tic the particle charge has not yet been determined 
from the Mie signal. Most particles of interest in astron- 
omy, astrophysics, atmospheric sciences, and laboratory 
experiments are however charged The particle 

charge is a rather important parameter. It determines 
the coupling of the particles among each other and to ex- 
ternal electro-magnetic fields. An optical measurement of 
it would be extremely useful. In principle, light scatter- 
ing contains information about excess electrons as their 
electrical conductivity modifies either the boundary con- 
dition for electro mag netic fields or the polarizability of 
the material 0, [9l4llj But how strong and in what spec- 
tral range the particle charge reveals itself by distorting 
the Mie signal of the neutral particle is an unsettled issue. 

In this Letter, we revisit Mie scattering by a nega- 
tively charged dielectric particle. Where electrons are 
trapped on the particle depends on the electron affinity 
X of the dielectric, that is, the offset of the conduction 
band minimum to the potential in front of the surface. 
For x < 0, as it is the case for MgO, CaO or LiF [I3 . ll3 j . 
the conduction band lies above the potential outside the 
grain and electrons are trapped in the image potential in- 
duced by a surface mode associated with the transverse 
optical (TO) phonon 14, [Hj]. The conductivity a s of 
this two-dimensional electron gas is limited by the resid- 
ual scattering with the surface mode and modifies the 
boundary condition for the electromagnetic fields at the 
surface of the grain. For x > , as it is the case for 
AI2O3 or SiC>2, electrons accumulate in the conduction 
band forming a space charge [Hj]. Its width, limited by 
the screening in the dielectric, is typically larger than a 
micron. For micron sized particles we can thus assume a 
homogeneous distribution of the excess electrons in the 
bulk. The effect on light scattering is now encoded in 
the bulk conductivity of the excess electrons ah, which 



is limited by scattering with a longitudinal optical (LO) 
bulk phonon and gives rise to an additional polarizabil- 
ity per volume a — -iTriab/uJ, where uj is the frequency 
of the light. We focus on the scattering of light in the 
vicinity of anomalous optical resonances which have been 
identified for metal particles by Tribelsky et al. 3, 17j |. 
These resonances occur at frequencies w where the com- 
plex dielectric function e(w) = e'(ui) + ie"(v) has e' < 
and e" -^i 1. For a dielectric they are induced by the 
TO phonon and lie in the infrared. Using Mie theory, 
we show that for submicron-sized particles the extinction 
resonance shifts with the particle charge and can thus be 
used to determine the particle charge. 

In the framework of Mie theory, the scattering and 
transmission coefficients connecting incident (i), reflected 
(r) , and transmitted (t) partial waves are determined by 
the boundary conditions for the electric and magnetic 
fields at the surface of the particle 0, [H[ . For a charged 
particle with x < the surface charges may sustain a 
surface current K which enters the boundary condition 
for the magnetic field. Thus, e r x (Hj + H r — Hj) = 
^f-K, where c is the speed of light 11]. The surface 
current K = <7 s Eii is induced by the component of the 
electric field parallel to the surface and is proportional to 
the surface conductivity a s . For x > the bulk surplus 
charge enters the refractive index N = y/e + a through 
its polarizability. Matching the fields at the boundary of 
a dielectric sphere with radius a gives, following Bohren 
and Hunt [ll], the scattering coefficients 
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where for x < (x > 0) the dimensionless surface con- 
ductivity t(w) = 47rer s (u;)/c (r = 0) and the refrac- 
tive index N — yfe (N — y/e + a), the size parameter 
p = ka = 27ra/A, where k is the wavenumber, ^p n {p) = 

V / 7rp/2J„ +1/2 (p), UP) = V 7T P/ 2H nll/2(P) With J n{p) 

the Bessel and Hn\p) the Hankel function of the first 
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kind. As for uncharged particles the extinction efficiency 
becomes Q t = -(2/p 2 ) J2n=i ( 2 « + 1) ^ e « + K)- Any 
effect of the surplus electrons on the scattering of light, 
encoded in a n and W n , is due to the surface conductivity 
(x < 0) or the bulk conductivity (x > 0). 

For x < we describe the surface electron film in a 
planar model to be justified below. For the dielectrics we 
consider, the low-frequency dielectric function is domi- 
nated by an optically active TO phonon with frequency 
ujto- For the modelling of the surface electrons it suffices 
to approximated it by e(oj) = l+u> 2 ^ (eo~ 1)/(lu^ — lo 2 ), 
where eo is the static dielectric constant, and allows 
for a TO surface mode whose frequency is given by 
e(uj s ) = -1 leading to uj s = uj TOy /(l + e )/2 [HI]. 
The coupling of the electron to this surface mode con- 
sists of a static and a dynamic part [20J. The for- 
mer leads to the image potential V = — k$e 2 / z with 
A = (e — l)/(4(e + 1)) supporting a series of bound 
Rydberg states whose wave functions read 0„k(x, z) — 
\J Ao /AaBTin\ 2 e ikx W ni i/2 (2A z/nas) with as the Bohr 
radius, k — (k x , k y ), x = (x, y), and A the surface area. 
Since trapped electrons are thermalized with the surface 
and the spacing between Rydberg states is large com- 
pared to ksT, they occupy only the lowest image band 
n = 1. Assuming a planar surface is justified provided 
the de Broglie wavelength XdB of the electron on the sur- 
face is smaller than the radius a of the sphere. For a 
surface electron with energy Ekm/ks = 300 A one finds 
XdB ~ 10~ 8 cto. Thus, for particle radii a > lOnm 
the plane-surface approximation is justified. The dy- 
namic interaction enables momentum relaxation paral- 
lel to the surface and limits the surface conductivity. 
Introducing annihilation operators Ck and oq for elec- 
trons and phonons, the Hamiltonian describing the dy- 
namic electron-phonon coupling in the lowest image band 



reads, H = J2w e kC k c k + Hlj s J2q a Q a Q + H m t [2l| with 
Hint = Z)k,Q(^k,Q/v / A)c k+Q (a Q - aL Q )c k , where the 
matrix element is given by (to is the electron mass) 
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Within the memory function approach [22j the surface 
conductivity can be written as 



e 2 n s 



m uj + M(uj) 



(3) 



with n s the surface electron density. Up to second order 
in the electron-phonon coupling the memory function 
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FIG. 1. Dimensionless surface conductivity r = r' + ir" for 
MgO for n B — 10 13 cm~ 2 (left) and polarizability of excess 
electrons a — a' + ia" for AI2O3 for rib — 3 x 10 17 cm~ 3 
(right) as a function of the inverse wavelength A -1 . 



where v = lo/lo s , S — /3huj s , 7 = y/2K^h/ a 2 B mui s , 
and I a {x) = r o °° dte-^^ t+t ^a e /(a + y/i) 6 which for low 
temperature, that is x — > 00, has the asymptotic form 
I a (x) ~ y/ 7T / 'xe~ 2x a 6 1 '(1 + a) 6 . Since M(ui) is indepen- 
dent of n s the surface conductivity is proportional to n s . 

For x > the bulk conductivity is limited by a 
longitudinal optical (LO) phonon with frequency ujlo- 
The coupling of the electron to this mode is described 

by H mt = Ek, q M 4+ q c k (a, + a f _ q ) /Wq Q , where 



M = y 2ire 2 foujLo (eoo 1 ~ e o X )- Employing the memory 
function approach, the bulk conductivity is given by Eq. 
(|3]) where n s is replaced by the bulk electron density 
rib and m by the conduction band effective mass to*, 
the prefactor of the memory function (Eq. (J3J)) is then 
4e 2 Vm*Lu L od(e^ ~ e^/^y/J^W), and 



= ^-|^ + l| e -^+ 1 )/ 2 A 1 (r5| J , + l|/2) 
e° — 1 

+ -^-rl^ " Ae-'^-WKi (8\u - l|/2) , (6) 
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where v — lo/lolo, 8 — flfiuJLO, and K\(x) is a modified 
Bessel function. For low temperature, i.e. 8 — > 00 ~ 

y/V/s^7TT\0(- v -i). 

To exemplify light scattering by a charged dielectric 
particle we will consider a MgO (x < 0) and an AI2O3 
(x > 0) particle 24j. The charge effect on scattering 
is controlled by the dimensionless surface conductivity 
t = t' + ir" (for x < 0) or the excess electron polariz- 
ability a — a' + ia" (for X > 0), both shown in Fig. [T] 
which are small even for a highly charged particle with 
n s = 10 13 cto~ 2 (corresponding to iij, = 3 x 10 17 cto -3 
for x > and a — l/im). The electron-phonon coupling 
reduces r" and a' compared to a free electron gas where 
M(lo) = 0, implying r' = and a" = 0. For T = OA, 
t' = (a" = 0) for A -1 < Xj 1 = 909cto-\ the in- 
verse wavelength of the surface phonon ( A -1 < X~[o = 
807cto _1 , the inverse wavelength of the bulk LO phonon) 
since light absorption is only possible above the surface 
(bulk LO) phonon frequency. At room temperature t" 
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FIG. 2. Dielectric constant e = e' + ie", refractive index TV = 
n + ik (top) and extinction efficiency Q t (bottom) depending 
on the particle radius a for MgO and AI2O3 as a function of 
the inverse wavelength A -1 . 
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FIG. 3. Magnification of the extinction resonances depending 
on A -1 and a. Crosses indicate their maxima for n s = (red), 
2 x 10 13 (green), and 5 x 10 13 cm -2 (blue) at T = 300K for 
MgO (left panel) and for nt, = 3n a /a for AI2O3 (right panel). 
Open circles indicate the maxima for free electrons. 



and a' still outweigh r 1 and a". The temperature ef- 
fect on t" is less apparent for A -1 > 300cm -1 than for 
a' but for A -1 < 300cm -1 a higher temperature lowers 
t" considerably. The upper panel of Fig. [5] shows the 
complex dielectric constant e and the refractive index N. 
For MgO we use a two-oscillator fit for e 25|, 26j ], In the 
infrared, e is dominated by a TO phonon at 401cm -1 . 
The second phonon at 640cm -1 is much weaker, justi- 
fying our model for the image potential based on one 
dominant phonon. Far above the highest TO phonon, 
that is, for A -1 > 800cm -1 (A -1 > 900cm -1 ) for MgO 
(AI2O3) e' > and e" <C 1. For these wavenumbers a 
micron sized grain would give rise to a typical Mie plot 
exhibiting interference and ripples due to the functional 
form of a r n and b r n and not due to the dielectric constant. 
Surplus electrons would not alter the extinction in this 
region because |e| ^> |t| and |e| 3> \a\. 

To observe a stronger dependence of extinction on the 
parameters e and r or a, we turn to 400cm -1 < A -1 < 
700cm -1 for MgO (700cm -1 < A -1 < 900cm -1 for 
AI2O3) where e' < and e" <C 1 allowing for optical 
resonances, sensitive to even small changes in e. They 
correspond to resonant excitation of transverse surface 
modes of the sphere 27J. For a metal particle the res- 
onances are due to plasmons and lie in the ultraviolet 
a dielectric the TO phonon induces them. 
As the polarizability of excess electrons, encoded in r or 
a, is larger at low frequency, the resonances of a dielectric 
particle, lying in the infrared, should be more suscepti- 
ble to surface charges. The lower panel of Fig. [5] shows a 
clearly distinguishable series of resonances in the extinc- 
tion efficiency. The effect of negative excess charges is 
shown by the crosses in Fig. [3] The extinction maxima 
are shifted to higher A -1 for both surface and bulk excess 
electrons. For comparison the circles show the shift for 
a free electron gas. The effect is strongest for the first 



resonance, where a surface electron density 10 cm (or 
an equivalent bulk charge) , realized for instance in dusty 
plasmas [28| . yields a shift of a few wavenumbers. 

The shift can be more clearly seen in Fig. 2] where 
the tail of the first resonance is plotted for MgO on an 
enlarged scale. The main panel shows the extinction ef- 
ficiency for n s = 10 13 cm -2 with its maxima indicated 
by blue dots. Without surface charge the resonance is at 
A -1 = 606cm -1 for a < 0.25pm. For a charged particle 
it moves to higher A -1 and this effect becomes stronger 
the smaller the particle is. The line shape of the extinc- 
tion resonance for fixed particle size is shown in the top 
and bottom panels for a = 0.2/im and a — 0.05pm, re- 
spectively. For comparison data for other surface charge 
densities are also shown. Figure 2] also suggests that the 
resonance shift is even more significant for particles with 
radius a < 0.01/im where the planar model for the image 
states is inapplicable. An extension of our model, guided 
by the study of multielectron bubbles in helium [29j j . 
requires surface phonons, image potential and electron- 
phonon coupling for a sphere. Due to its insensitivity to 
the location of the excess electrons, we expect qualita- 
tively the same resonance shift for very small grains. 

As we are considering particles small compared to A we 
expand the scattering coefficients for small p. To ensure 
that in the limit of an uncharged grain, that is, for r — > 0, 
a r n and b r n converge to their small p expansions [18| . we 
substitute t = r/p before expanding the coefficients. Up 
to 0(p 3 ) this yields a{ — a r 2 — b r 2 — and only b\ ~ 
0(p 3 ) contributes. Then the extinction efficiency, 



Qt 



12p(e" + a" + 2t'/ P ) 



(e' + a' + 2 - 2t"/ P ) 2 + (e" + a" + 2r'/p) 2 



, (7) 



where the excess charges enter either through r with a — 
for x < or through a with r = for \ > 0. For 
r, a —> this gives the limit of Rayleigh scattering. The 
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FIG. 4. Middle panel: Extinction efficiency Qt as a function 
of the inverse wavelength A" 1 and the radius a for a MgO 
particle with n s = 10 13 cm~ 2 and T = 300K. The dotted 
lines indicate the extinction maximum for n s — (black), 
2 x 10 12 (green), 5 x 10 12 (red), and 10 13 cm" 2 (blue) obtained 
from JSJ. Top and bottom panel: Extinction efficiency Qt 
as a function of A -1 for a — 0.2/im (top) and a = 0.05fim 
(bottom) for different surface electron densities. 



resonance is located at the wavenumber where 



2 - 2r"/p = 



(8) 



and has a Lorentzian shape, already apparent from Fig. 
HI provided e" and r' (or a") vary only negligibly near 
the resonance wavelength. For an uncharged surface the 
resonance is at Aq 1 for which e' = —2. For x < the 
shift of the resonance is proportional to r" and thus to n s , 
provided e' is well approximated linearly in A -1 and t" 
does not vary significantly near Aq 1 . In this case, we sub- 
stitute in © the expansions e' = — 2 + c e (A~ 1 — A^ 1 ) and 
t" = c T n s where c e = g ^_i and c T = 2i_| A -i. Then 

the resonance is located at A -1 = A^ 1 +c T rt s /(7rc e aA ~ 1 ). 
For x > the resonance is located at A -1 = A^ 1 — c Q nf,/c e 



where c a = f-| 



The dotted lines in Fig. [4] give the 



location of the resonance obtained from Eq. © for MgO, 



where A 



606cm for several surface electron den- 
sities. They agree well with the underlying contour cal- 
culated from the exact Mie solution, as exemplified for 
The proportionality of the resonance 



= 10 13 cm 



shift to n s (rib) can also be seen in Fig. [S]where we plot on 
the abscissa the shift of the extinction resonance arising 
from the surface electron density given on the ordinate 
for LiF 0, MgO (x < 0) and A1 2 3 ( X > 0). Both bulk 
and surface electrons lead to a resonance shift. To illus- 
trate the similarity of the shift we consider ([5]) for free 
electrons, which then becomes e' — 2N e e 2 / (ma 3 u> 2 ) = —2 
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FIG. 5. Position of the extinction resonance depending on 
the surface charge n 3 for LiF , MgO and AI2O3 (for equivalent 
bulk charge rib = 3n 3 /a) particles with different radii a. Solid 
(dashed) lines are obtained from the Mie contour (Eq. ©). 



for x < and e' ~3N e e 2 /(m*a 3 uj 2 ) = -2 for X > 0; N e is 
the number of excess electrons. The effect of surface elec- 
trons is weaker by a factor 2m* /3m where the 2/3 arises 
from geometry as only the parallel component of the elec- 
tric field acts on the spherically confined electrons. Most 
important, however, r/p and a enter the resonance con- 
dition on the same footing showing that the shift is es- 
sentially an electron density effect on the polarizability 
of the grain. We therefore expect the shift to prevail also 
for electron distributions between the two limiting cases 
of a surface and a homogeneous bulk charge. 

To conclude, our results suggest that for dielectric 
particles showing anomalous optical resonances the ex- 
tinction maximum in the infrared can be used to deter- 
mine the particle charge (see Fig. [5]). For dusty plasmas 
this can be rather attractive because established methods 
for measuring the particle charge [30 32] require plasma 
parameters which are not precisely known whereas the 
charge measurement by Mie scattering does not. Parti- 
cles with surface (negative electron affinity x, e -g- MgO, 
LiF) as well as bulk excess electrons (x > e.g. AI2O3) 
show the effect and could serve as model systems for sub- 
micron sized dust in space, the laboratory, and the atmo- 
sphere. These particles could be also used as minimally 
invasive electric probes in a plasma, which collect elec- 
trons depending on the local plasma environment. De- 
termining their charge from Mie scatte ring and the forces 
acting on them by conventional means [30h32| | would pro- 
vide a way to extract plasma parameters locally. 
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